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Abstract 

We consider the supersymmetric seesaw mechanism induced by the exchange of heavy 
SU(2)w triplet states, rather than 'right-handed' neutrino singlets, to generate neu- 
trino masses. We show that in this scenario the neutrino flavour structure tested at 
low-energy in the atmospheric and solar neutrino experiments is directly inherited from 
the neutrino Yukawa couplings to the triplets. This allows us to predict the ratio of 
the r — > /U7 (or r — > e~f) and fi — ► ej decay rates in terms of the low-energy neu- 
trino parameters. Moreover, once the model is embedded in a grand unified model, 
quark-flavour violation can be linked to lepton-flavour violation. 



1 Introduction 



Nowadays we have one more important piece of information concerning the flavour structure 
of the Standard Model: flavour violation is also present in the lepton sector. So far this 
has only shown up in the atmospheric and solar neutrino experiments [|T], The anomalies 
observed in these experiments can be interpreted in terms of neutrino oscillations which are 
the result of non- vanishing neutrino masses and mixing angles ||, f|. In the framework of the 
Standard Model (or of its supersymmetrized version) the latter properties effectively arise 
from the following lepton- number (L) violating d — 5 operator |J: 

' Y i J(L i H 2 ) {LjH 2 ) , Y V = Y T U (1) 
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where i,j = 1,2,3 are family indices, Ml is the energy scale where L is broken, Li are the 
SU{2)w lepton doublets and H 2 is the Higgs doublet with hypercharge Y = 1/2. Upon 
breaking SU{2) W x t/(l)y by the vacuum expectation value of the Higgs field, (H 2 ) = v 2 = 
vsmfi (v = 174GeV), the operator ([[]) induces Majorana masses for the neutrinos, 

v 2 

m ij _ J2_Y ij (2) 

Therefore, from here we can easily understand the origin of the tiny neutrino masses, as 
they may be suppressed by some large scale Ml- Taking for instance a neutrino mass 
m v ~ (10 _1 — 10~ 2 ) eV as indicated by the atmospheric neutrino data (and assuming a 
hierarchical neutrino spectrum) ||, the magnitude of Ml can be inferred, 

Y~ X M L ~ 10 15 GeV. (3) 

In the basis in which the charged-lepton Yukawa matrix Y e is diagonal, all the lepton-flavour 
violation is contained in the coupling matrix Y„, i.e. in the neutrino mass matrix: 

m u = "iTm^U 1 ", = diag(mi,m 2 ,m 3 ) (4) 

where mi, m 2 , are the neutrino mass eigenvalues, and the unitary matrix U is the lepton 
mixing matrix that appears in the charged lepton current £^(1 — 75)U v and is responsible 
for neutrino oscillations. 

In principle lepton-flavour violation could also be tested in other processes, such as — ► 
ej, t — > /i7 and r — > cy. The present experimental limits on these decays are M: 



BR(fi - 


-> erf) 
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1.2 x 10 


BR(t - 


■> ni) 
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1.1 x 10 


BR(t - 


-> ery) 
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2.7 x 10 



-11 



(5) 



which could be significantly improved, the first down to 10~ 14 and the second to 10 -9 
|§. In the Standard Model it is very hard to obtain interesting results because the decay 
amplitudes are strongly GIM-suppressed by the tiny neutrino masses, m 8 M w fIQ| . On 
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the other hand, in the framework of the minimal supersymmetric extension of the Standard 
Model (MSSM), those processes can be enhanced through the one-loop exchange of super- 



partners, if the masses of the latter are not too heavy and do not conserve flavour fll 
Concerning the latter property, we could expect that the underlying flavour theory dictates 
both the flavour structure of the standard fermion mass matrices and that of the correspond- 
ing supersymmetric scalar partners (see e.g. ||12|| ). Even if sfermion masses are universal (i.e. 
flavour-blind) at high energy, as in minimal supergravity or gauge mediated supersymmetry 
breaking scenarios, nevertheless flavour conservation can be broken in the sfermion masses 



by radiative effects due to flavour- violating Yukawa couplings |13], [14]]. In particular, the 
interactions that generate the operator (H) also induce lepton-flavour violation (LFV) in the 
slepton mass matrices by renormalization effects. A well-known and investigated example is 



that of the standard seesaw mechanism |L5| in which LFV is induced radiatively by Yukawa 



couplings of the SU(2) ^-doublet neutrinos with singlet neutrinos N (often referred to as 



right-handed neutrinos) [|T^, |16 |. Here, we would like to discuss another example of seesaw 
scenario in which the d = 5 operator (|I|) is obtained by exchanging heavy SU(2)w triplets 
with non-zero hypercharge. We recall that models with scalar triplets to generate Majo- 
rana neutrino masses have been considered in the literature for a long time, though they 
have received less attention than the standard seesaw scenario. For example, a model with 
spontaneous L-breaking was proposed in |l7j and later on extended |18| |. A triplet-exchange 



seesaw realization with explicit L-breaking was also introduced [19[], in a non-supersymmetric 
frameworkF]. A supersymmetric version of the latter scenario was recently introduced to have 
baryogenesis through leptogenesis |22| . In this work we shall further elaborate the supersym- 
metric triplet seesaw scenario. In particular, we shall discuss how LFV is radiatively induced 
in the slepton masses and show that this scenario is potentially more predictive than the 
N- induced seesaw one. 

This paper is organised as follows. In Sect. 2 we review both the standard see-saw 
mechanism and that in which triplet states are exchanged. This presentation will help in 
showing the differences between the two scenarios and outline the one-to-one correspondence 
between the neutrino parameters and LFV in the soft-breaking parameters that characterizes 
the triplet seesaw scenario. In Sect. 3 the triplet seesaw is embedded in SU(5) context. In 
Sect. 4 we show the general pattern of Y„ as derived from the low-energy neutrino data. 
In Sect. 5 we discuss the flavour- violation induced in the soft-breaking terms by radiative 
effects in the energy range above the triplet mass threshold. The renormalization group 
equations (RGEs) relevant to our study are confined in Appendix. In Sect. 6, after giving 
the qualitative behaviour of the £j — > iff branching ratios, we also discuss some numerical 
examples. We give our conclusions in Sect. 7. 



1 Recently, this scenario was studied for leptogenesis [|20| . An alternative seesaw mechanism, obtained by 
exchanging heavy SU(2)w triplets with zero hypercharge, was discussed in fl2l| . 
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2 Singlet versus Triplet seesaw 



First, we briefly review the standard seesaw mechanism in which singlet states N are ex- 
changed [[15] . The relevant superpotential terms at the scale where lepton number is broken 
are: 

Y%N.L,H 2 + ^N,N„ (6) 

where i,j = 1, 2, 3 are family indices, Y^ is an arbitrary 3x3 matrix of Dirac-like Yukawa 
couplings, while M^v is a 3x 3 symmetric mass matrix describing Majorana masses for the 
singlets N. If the N states are assigned lepton number L = — 1, the second term in (|6]) 
explicitly breaks L. After decoupling the heavy states N at the (overall) scale Ml, the 
lepton-number violating d = 5 operator (|1|) is generated, where Mj}Y v is identified as 
follows: 

Finally, at the electroweak scale Majorana masses for the neutrinos are obtained, 

m£ = ^Y? = v\ Y™]Vtf w Y#. (8) 

As a matter of fact, there is a certain degree of ambiguity in the effective neutrino mass 
matrix: its flavour structure reflects both that of the arbitrary Dirac-like matrix Y ^ and 
that of the matrix M^v- From a bottom-up perspective, this implies that neutrino masses 
and mixing angles, inferred by the low-energy neutrino data ||, may only reflect the effective 
Yukawa matrix Y u and the overall mass scale Ml (modulo radiative corrections) but cannot 
be unambiguously related to the more fundamental quantities, Y N and (for a recent 
discussion on this aspect see for example [^] and references therein). In other words, the 
low-energy parameters, described by Y v , which amount to 6 real parameters + 3 phases, are 
less than the number of the independent 'fundamental' physical parameters in Y^ and Mat, 
which instead are 12 real parameters + 6 phases. 

Now, let us consider the triplet seesaw scenario. The relevant superpotential terms are 

-^Y^LiTLj + -^=X X H X TH X + ^=X 2 H 2 TH 2 + M T TT + ^H X H 2 , (9) 

where the super-multiplets T,T are in a vector-like SU(2)w x U(l)y representation^, T ~ 
(3,1) and T ~ (3,-1), and H x is the Higgs doublet with hypercharge Y = —1/2. The 
matrix Y^ is in general a 3 x 3 symmetric matrix, Yj- = Y^\ If we assign lepton number 
L = —2 (2) to the triplet T (T), we can see that the Ai, A2-couplings explicitly break L. If 
instead we assign L = — 2 to T and L = to T, then the L-breaking parameters are Mt and 
Ai. In eq. @ the triplets T, T are represented as 2 x 2 matrices, namely 

/ T° -T + \ ( f — Lj 1 - \ 

T = (ia 2 )T.<r = ^ _^ T+ _^ ++ J , T = (ta 2 )T-a = ^ __^ f _ f f0 j , (10) 



2 Notice that in the supersymmetric picture two triplets, T and T, are required. 
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where the matrices a a (a — 1, 2, 3) are the Pauli matrices and the components Ti, T2, T3 and 
Ti,T2,T 3 of T and T, respectively, can be easily inferred?]. By decoupling the triplet states 
at the scale we again obtain the d = 5 effective operator of eq. dH). Now, however, the 
flavour structure of the matrix Y u is the same as that of Y^, that is: 



Mr 



X2 



- I rp , 



(11) 



and the Majorana neutrino mass matrix is given by: 



m 



Mr 



r '.l 



V2X2 

M T 



(12) 



Fig. [l| shows the diagrams inducing the supersymmetric d = 5 operator (|T]) in the iV-seesaw 
(a) and T-seesaw (b). In terms of component fields, the neutrino mass operator is generated 
through the exchange of the N fermion component in the former realization, while in the 
latter realization it is generated by the exchange of the T scalar component and of the T 
F-component which gives rise to the effective scalar coupling \2MtH2T<H2- This latter 
feature also allows us to appreciate the supersymmetric version @ of the T-induced seesaw. 
In the non-supersymmetric case [pi the coefficient of the cubic interaction HT^H is an 



independent mass parameter, say A, and therefore m p 



Ho 



Ho 



N 



(b) 



T 
T 



Figure 1: Contributions to the LLH2H2 effective operator: (a) heavy singlet exchange; (b) heavy 
triplet exchange. 

The difference with the standard iV-induced seesaw case is manifest: in the T-seesaw 
scenario the neutrino mass matrix can be directly related to the fundamental matrix and 
to one relevant mass scale parameter, Mt/\2- I n particular, the amount of lepton-flavour 
violation measured at low-energy through the lepton mixing can directly be linked to the 
only source of lepton-flavour violation, the Yukawa matrix given at the scale Mt- The 
counting of the independent parameters reveals indeed that in this model the independent 
Yt parameters, 6 real parameters + 3 phases, are just matched by the low-energy physical 
parameters. There are left other 2 real parameters, Mt and A2, whose ratio is directly 
involved in the neutrino mass generation, plus Ai which in general is complex and is not 
directly connected to the neutrino parameters. In summary: while in the iV-induced seesaw 

3 As a consequence of the representation adopted for T and T ([H]), the 5 l [/(2)vi/-invariant mass term in 
cq. (§) is to be understood as M T TT = M T Ti(Tia 2 fi(7 2 ). 
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scenario there are two lepton- flavour violation sources, the matrices and Mjv, in the 
T-induced seesaw scenario the matrix is the only source of lepton-flavour violation^. 

The implications of all this are even more dramatic when we consider the lepton-flavour 
violation induced through renormalization effects. Since a long it has been pointed out 
that the Yukawa couplings Y^v in eq. (|6]) can induce non-vanishing lepton-flavour violating 
entries in the mass matrices of the left-handed sleptons through radiative corrections | I3fl , 



even in the minimal SUSY scenario with universal soft-breaking terms at the GUT scale 



Mi 



G, m\ 



moll. The form of the LFV entries is 



m|)y oc mo(Yj v -Y 7V )ylog 



M. 



G 



M 



N 



(13) 



In this scenario, according to our previous discussion, the size of LFV cannot be unambigu- 
ously predicted in a bottom-up approach making use of the low-energy dataf]. 

Now let us see what can occur in the T-seesaw scenario. In this case, the lepton-flavour 
violating entries are directly connected to the effective neutrino mass matrices, as 



[xa\)ij oc mo(Y^Y T )ylog 



Mi 



G 



M 7 



(14) 



or more explicitly, 



M T 
\ 2 v% 



Mg 
M T 



M T 

A 2 t>2 



U(m^ ) ) 2 U t 



log . 

ij M T 



(15) 



This expression enables us to univocally predict the ratio of the lepton-flavour violation in 
the 2-3 sector with that in the 1-2 sector, essentially in terms of the low-energy parameters, 
namely: 



U(m^) 2 Ut 



m 



l)He 



or for the 1-3 sector 



[U(m^UV 



m 



Thus we can relate the rate of the r — > /ry or r 
This is the main feature of our discussion. 



e7 decay with that of the \i 



(16) 

(17) 
e7 decay. 



4 Notice that in this seesaw scenario one triplet pair T, T is enough to generate non- vanishing mass for 
all three neutrinos, whilst in the standard seesaw three singlets N are necessary for that. Also notice that, 
from the flavour point of view, the seesaw realized by exchanging hyperchargeless triplets T' ~ (3, 0) [ pl[ is 
more similar to the iV-induced seesaw. Indeed: i) such triplets T" are exchanged in the same channel as the 
singlets N (see Fig. [l] (a)); ii) three hyperchargeless triplets are required to give mass to all three neutrinos; 
iii) two sources of flavour-dependence appear, the matrix Yt* in the Yukawa interaction ~YtLT' H2, and the 
mass matrix Mp of the triplets. 

5 The authors of |2j| indeed regard the combination YjyYjv, responsible for the LFV in the slepton 
masses, as a further 'observable' which provides the lacking 6 real parameters + 3 phases necessary to fully 
determine both Yn and Mjv. 
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This scenario is susceptible of being further elaborated. Indeed, the presence of the extra 
SU( 2) w -triplet states T, T at intermediate energy would spoil the gauge coupling unification 
which takes place with the field content of the MSSM. A simple way to save gauge coupling 
unification is to introduce more states X, to complete a certain representation R - such that 
R = T + X - of some unifying gauge group G, G D SU (3) x SU(2) W x U (1) Y ■ Thus we can 
envisage three (minimal) scenarios: 

(A) - The states X, though with mass M x ~ M T , are assumed to have vanishing or negligible 

interactions with the other states, except for the gauge interactions. This may be the 
case either if we prefer not to embed the theory in any definite grand unified theory 
(GUT) or if we embed it in some GUT and the Yukawa interactions of the states T and 
X have different strength due to GUT breaking effects. In particular, we could have 
negligible Yukawa coupling for the fragments X and, on the contrary, non-vanishing 
Yukawa coupling for the fragments T. 

(B) - The theory is embedded in a GUT but, contrary to the ansatz (A), the Yukawa 

couplings of the states X are assumed to be non-vanishing and related to those of the 
triplet partners T. Indeed, this is generally the case in minimal GUT models. In this 
case we will generate not only lepton- flavour violation but also closely related flavour 
violation in the quark sector (related to the X-couplings). 

(C) - There are no extra states X or, equivalently, they are considered to be split in mass 

from the triplets T and be decoupled at a scale fi > Mq- In this case, the simple 
unification of gauge couplings is lost and large threshold corrections are needed to 
recover it. 

3 SUSY SU(b) scenario with SU(2) W triplets 

As already mentioned in the previous section, the extra states T, T with mass Mt much 
below the GUT scale would destroy the gauge coupling unification. In principle the latter 
property could be recovered at the price of invoking large threshold corrections. In the next, 
however, we prefer to maintain the simple gauge coupling unification. To this purpose, the 
field content of the model can be minimally extended by adding the other components of 
the ££7(5) representations, 15 and 15, in which the triplets T and T can indeed fit. In terms 
of £77(3) x SU(2) W x U(l)y representations, the 15-multiplet decomposes as follows: 

15 = £ + T + Z, 

£~(6,l,-jj), T~ (1,3,1), Z~(3,2,i), (18) 

(the 15-decomposition is obvious). The presence of these extra states fitting a complete 
GUT multiplet changes the value of the gauge coupling a G at the GUT scale, with respect 
to the MSSM case, but does not modify the value of the unification scale M G (to one loop 
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accuracy). The /3-functions of the gauge couplings in the RGEs get modified as follows 
(a = 1,2,3): 



16vr — = B a9a , 

3 8 1 7 

Bi = b x + -(-n s + 3n T + -n z ) = h + -n 15 , 
5 3 o 2 

3 7 
B 2 = b 2 + 2n T + -n z = b 2 + -n 15 , 

5 7 
B 3 = b 3 + -ns + n z = h + ^n 15 , (19) 

where b a are the coefficients of the /3-functions in the MSSM, namely b\ — y , b 2 — 1, 63 = —3, 
and we have explicitly shown the contribution of the new states (ns — Ns + N§, = 
N T + Nf, n z = N z + N z and n 15 = N 15 + Nf 5 in a self-explanatory notation). As expected, 
for each B a the overall contribution of S, T, Z just reproduces the Dynkin index | of the 
SU(5) 15 representation. The enhancement of the /5-functions makes the gauge couplings 
increase faster at higher energy. For instance, by using the low-energy values of a em , a s and 
sin 2 Ow and assuming an average SUSY threshold close to the top mass, for T, S, Z masses 
around 10 14 GeV we find that at one-loop g\ and g 2 get unified at Mq ~ 2 x 10 16 GeV to the 
common value go ~ 0.88, while the value of #3 differs by one per mill or so. In the MSSM 
we would find go ~ 0.71. 

The SU(5) invariant superpotential (omitting for simplicity the flavour indices) reads as: 

W sm = "^YisS 15 5 + -^Ax 5 H 15 5 H + X 2 5 h 15 5 H + Y 5 10 5 5^ 

+Y 10 10 10 5 H + M 15 15 15 + M 5 5 (20) 

where the matter multiplets are understood as 5 = (d c ,L), 10 = (u c ,e c ,Q) and the Higgs 
doublets fit with their coloured partners, t, £ as 5# = (t, H 2 ), 5# = (i, Hi). In the 15-multiplet 
the states S, T and Z are accommodated as: 

(gab _2_Z a j \ 
— Z U T ij ' J 

where the S'?7(5)-indices A, B = 1, 2, 3, 4, 5 are decomposed into SU (3) indices a,b — 1,2,3 
and SU{2)w indices i,j = 4, 5, A = (a,i),B = (b,j). In eq. fl2jj ) it is understood that 

S aa = £aa ? gab = J_gab ( Q ^ Q ^ rpii = fii^ rpj = _^fij ( • ^ ( cfr gq ^ where the 

fields S ab , T l i are those canonically normalized. It is well-known that the minimal SU(5) 
model in which the Yukawa matrices Y 5 , Y 10 are true constants is not phenomenologically 
satisfactory. The latter should be rather understood as field dependent quantities, e.g. 

Y 5 ($) = Yj 0) + Y^ 1} $/M H where $ is the adjoint 24 of SU(5) and M is some cutoff 

scale larger than the GUT scale Mq. This perspective allows us to correct certain SU(5)- 
symmetry relations, such as Y^ = Yj pi. Moreover, some mechanism is also necessary to 
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split the masses^ of the doublet Hi j2 and triplet t, t components of 5#, 5h not to have fast 
proton decay mediated by the coloured states t, t |25j. We shall therefore adopt this point of 
view for the whole SU(5) extended model of eq. (|20|) . In the 5'f/(5)-broken phase, beneath 
M G , the superpotential reads as: 

-^={Y T LTL + Y s d c Sd c ) + Y z d c ZL + Y d d c QH 1 + Y e e c LH 1 + Y u u c QH 2 
v2 

+^=(X 1 H 1 TH 1 + \ 2 H 2 TH 2 ) + M T TT + M Z ZZ + M S SS + ^iH x H 2 . (22) 
v2 

As already mentioned, the couplings involving the coloured triplets t, t do not appear as the 
latter are assumed to decouple at the GUT scale Mq to suppress dangerous B — L violating 
d = 5 operators. On the contrary, at the decoupling of the 15 fragments, no B — L violating 
d = 5 operators are induced, apart from the 'neutrino' operator, since the 15 states do not 
couple to the matter multiplets 5, 10. Only flavour- conserving d = 6 operators are generated, 
i.e. (Ld c )(Ld c ), (d c d c )(d c d c ), (LL) (LL) in the Kahler potential, which, being suppressed by 
the square of the large scale M15, are not relevant for the low-energy phenomenology. Notice 
that in the minimal case the masses Mr, Ms, Mz are equal at the GUT scale, 

M s = M T = M z = M° 5 . (23) 

However, this unification relation could be modified e.g. due to ^-insertions, as mentioned 
abovefj. Similarly, the unification of the Yukawa couplings at the GUT scale, 

Y S = Y T = Y Z = Y? 5 , (24) 

could either hold or not. Finally, we stress that in this SU(5) framework the flavour violation 
is encoded not only in Yy (as already elucidated above) but also in Y^ and Y z- Therefore, 
non- vanishing barion-flavour violating entries in the mass matrix ra?~ of the sdown squarks 

d c are induced by radiative corrections. In summary, the three scenarios (A), (B) and (C), 
put forward in Section 2, can be rephrased as follows: 



(A) - All fragments T, 5* and Z have the same mass, i.e. eq. ( |2"B| ) holds at Mq- However, 

the couplings Y^, Y z are assumed to be negligible, i.e. eq. (|24D does not hold. In 
this case only the interactions with the triplets T and so the couplings Y? drive the 
lepton-flavour violation in the slepton scalar masses. 

(B) - Both the masses and the Yukawa couplings of the 15 states are unified, i.e. both 

eqs. (^) and ( p4|) hold. Therefore all the couplings Y^, Yt, Y z will induce flavour 
violation in both the slepton L and squark d c masses. 



6 Also for the minimal 'technical' realization of the doublet-triplet splitting [g5| we can give an inter- 
pretation of the mass parameter M5, analogous to that adopted for the Yukawa couplings, i.e. Ms(<!>) = 
M 5 ° + A 5 $ + • • •. 

7 Moreover, even if eq. (|2^) holds at Mq, renormalization effects split the masses at lower energies. 
However, the relative splitting is not large and we will decouple all the components T, S, Z at the common 
threshold scale Mt- 
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(C)- The triplet mass Mt is much smaller than Ms and Mz, which are 0(Mq)- This 
could be achieved for instance by tuning the coefficients of the singlet and adjoint 
components of the 15 15 'mass', in analogy to what is done for the doublet-triplet 
splitting for 5#, 5#. 

In the following we shall focus on (A) and (B) and disregard the case (C). 



4 Yt from neutrino masses and mixing angles 

Now, we relate the low-energy parameters with the relevant neutrino Yukawa couplings by 
adopting a bottom-up criterion. By decoupling the states T, T, the d — 5 effective operator 
emerges 

^-Y^(L^ 2 )(L,if 2 ), (25) 

where the matrix Y^, through the matching of eq. (|TT|) , can be connected to Y u which 
parameterizes the usual d = 5 operator ([I]) Q. 

We recall that the data from solar and atmospheric neutrino experiments concern the 
neutrino mass eigenvalues mi,m 2 ,rri3 and mixing angles. Therefore, in the basis in which 
the Yukawa matrix Y e is diagonal, eqs. (§) and (f|) allow us to determine the coupling matrix 
Y u /Ml at low-energy and then, taking into account the running up to Mr PB]| , also at the 
scale The mixing matrix U is parameterised in the standard way: 

C12C13 S12C13 Si3e~ J<5 

U = j -S12C23 - c 12 s 2 3Sne l5 ci 2 c 23 - Sns^s^e 15 S23C13 ! (26) 
S12S23 ~ c 12 c 23 s 13 e lS -C12S23 - s 12 c 2 3Si3e lS c 23 ci 3 

where and c^- are the cosine and sine respectively of the three mixing angles 6i 2 , 6 2 3, 613 
and S is the CP-violating phase which in the following is neglected for simplicity. As for 
the phenomenological input, we assume maximal 2-3 mixing, #23 = 45°, as required by the 
atmospheric neutrino data 0. We consider hierarchical neutrino mass spectrum, mi <C 
m 2 <C 777,3, sucri that it reads as: 



mi 



0, m 2 = (AmL) 1/2 , m 3 = (Am 2 tm ) 1/2 , (27) 



where we take Am 2 tm ~ 3 x 10~ 3 eV 2 . As regards the solar neutrino case, the most favoured 
range for Am^ is that selected by the large-mixing angle (LMA) solution, Am^ ~ 6 x 
lCT 5 eV 2 ||. The corresponding best fit value for the mixing angle is 6*12 ~ 33°. However, for 
the sake of discussion in the following we bear in mind also different possibilities, such as the 

8 The expression given in eq. ( p5| ) is the leading contribution to the neutrino mass operator. This arises 
from the T F-term scalar interaction, X^MtT* The H\ F-term scalar interaction, /iXiH^THi, gives 
rise to the sub-leading contribution \\([ji/Mt){v2Vx)/Mt- 

9 The radiative corrections from the electroweak scale up to Mt are not important in the case of hierar- 
chical neutrino mass spectrum and amount to an overall common factor ^7j. Nevertheless these effects are 
incorporated in the numerical analysis. 
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typical values Am s 2 ] ~ 5 x 10~ 6 eV 2 and 9 12 ~ 4° of the small-mixing angle (SMA) solution. 
Taking also into account the CHOOZ limit, sin 6*13 < 0.1 |[28|| , we set #13 = for simplicity 
and later comment also on the non-zero #13 case. Then at low-energy the symmetric matrix 
Y u appears as: 



Y M ' 



v 2 



( m 2s\ 2 ^m 2 c 12 s 12 ~^m 2 c 12 s 12 

i(m 3 + m 2 c\ 2 ) |(m 3 -m 2 cf 2 ) |. (2? 
V |(m 3 + m 2 c 2 2 ) 

By considering the phenomenological inputs, we observe that the entries of the 2-3 sector 
are all comparable and can be of order one (essentially irrespectively of the 9 12 value) if the 
scale Mi is close to Mq ~ 10 16 GeV. The remaining entries are one order of magnitude 
smaller than those of the 2-3 sector, for 9\ 2 in the LMA range, while they are much smaller 
for 9 12 in the SMA range. Up to an overall factor due to the renormalization effect, the 
structure obtained for Y u is finally transferred to Y T at the scale M T according to eq. (|TT|). 
We stress again that in this T-seesaw scenario the bottom-up approach here adopted is the 
most general one since the structure of Yy is unambiguously fixed by the experimental data 
themselves. What can be left to our choice is the overall scale Mt/X 2 - We shall take Mt in 
the range 10 11 -j- 10 15 GeV and vary A2 in an appropriate range. 

5 Lepton-flavour violation in the soft-breaking terms 

The general soft SUSY-breaking terms in our model are given as 10 : 

- £ soft = Dm 2 L L + e c m|e ct + d c m%d c1 + m 2 H Il\R x + m 2 H2 H\H 2 
+ (#ie c A e L + Hid c A d Q + ^M a A a A a + BiiH x H 2 + h.c.) 



+ 



-^={A T LTL + A s d c S<i c ) + A z d c ZL + -^=(A 1 H 1 THi + A2H2TH2) 



+B T M T TT + B S M S SS + B Z M Z ZZ + h.c. 

Z^Z + m\Z^Z (29) 

where we have shown in the first lines the relevant MSSM terms, according to the standard 
notation (the soft mass terms for the sleptons, squarks, Higgs bosons, the trilinear terms 
and the gaugino masses M a ), while in last lines we have collected the new terms involving 
T, S, Z. In the following we assume at the GUT scale, irrespectively of the scenario (A) or 
(B)0 

2 2 2 2 11 
mg = m e - c = = m Q Jl, 



10 For the sake of simplicity, in the following we disregard all what concerns the up-quark sector parameters, 
such as the scalar masses mj~ etc. because they do not enter directly our discussion. 

11 These universal boundary conditions are mainly motivated by simplicity and are not assumed for other 
soft parameters. For instance, we do not declare the soft-breaking mass of the doublet scalar H2, as it does 
not directly enter in our analysis. In this way the /x-parameter in the superpotential is not constrained by 
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m s = m s = m T = nif = m z = m z = % = m , 

A e = A Y e , A d = A Y d , A l = A X U A 2 = A X 2 , 

M. M-2 M 3 M g , (30) 

where mo is the universal scalar mass, Aq is the universal mass parameter for the trilinear 
terms and M g is the common gaugino mass. As for the remaining trilinear couplings, their 
GUT conditions are in the scenario (A): 



and in the scenario (B) 



A Y T , A S = A Z = 0, (31) 



A r = A s = A z = A Y T . (32) 



The matrix Yt which appears in eqs. fl3T| - |3"2"D is understood to be evaluated at the GUT 
scale. Once Yt is determined at the scale Mt by the low-energy data as we have seen above, 
its evolution from Mt up to Mq is given by: 

16vr 2 ^ = Y r (-^-7 S | + YtY e + 6Y|,Y T + '&(Y r Yt) + 3YiY z + |A 1 | 

+ (YjY: + 3Y|Y^)Y T . (33) 

(See also Appendix.) In the scenario (B) the Y T RGE is also coupled to those of Yg, Y z for 
which the initial conditions at the scale Mt are determined iteratively under the constraints 
of eq. fl24|). Below Mq the universal pattern (|30|) of m|,m 2 - c etc. is spoiled by radiative 



effects induced by Yt, Y§, Y z . Then we have to evaluate the soft-breaking parameters at 
low-energy by solving the corresponding RGEs. These have been computed at one-loop and 
collected in Appendix. In the leptonic sector we need to know the SUSY breaking matrices 
m|,m| c and A e to finally compute the LFV decay rates. They all receive flavour blind 
corrections from the gauge interactions which do not alter the flavour-conserving structure 
they have at the GUT scale (see eq. ((50)). However, they can acquire LFV entries (i.e. off- 
diagonal entries) if they get radiative corrections from the LFV Yukawa-matrices Yt, Y z . 
In the leading-log approximation, and neglecting radiative corrections induced by Y e , Y^, in 
the picture (A) the LFV entries at low-energy are given by (i ^ j): 



« ^(9m^ + 3^)(Y^YT) u log MG 

2 



TT 2 ' U U " 1 ' J ° Mn 



(Ae)ii « ^i4o(Y e Y},Y r )ylog^. (34) 



the electroweak radiative-breaking condition and will be fixed independently. We would like to make clear, 
however, that the aim of this work is to present the global features of the SUSY T-induced seesaw. Therefore, 
the choice of these initial conditions, as well as of other parameters such as tan/3 and /i, which enter in the 
computations of the decay rates, is only made for illustrative purposes. 
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and in the picture (B) we find: 







-1 


(ml 


1 y 


8tt 2 




y 


« 0, 






-9 


(a.: 


y 


8tt 2 



;i8m2 + 6^)(Y^Y T )^log^, 



+ Mr- 

A (Y e Y^Y T )ylog^, (35) 



and in the squark sector: 



n !)y « g^(18m2 + 6^)(Y}.Yr)«log^, 

A„)« « f^A)(Yj,Y r Y d )„log^, (36) 



where we have taken into account the S77(5)-universality for the Yukawa matrices (24). So 



the mass matrix m~ c remains diagonal and then flavour- conserving, while both m? and A e 
acquire LFV elements. Once the low-energy neutrino observables are fixed, the magnitude 
of these LFV elements will depend on the matrix Y^Y-p ~ (Mr/A2f|) 2 mJ,m i/ , that is on the 
triplet mass threshold and on the coupling constant A2. Thus the relative size of LFV in 
the 2-3 sector and 1-2 sector can be approximately predicted in terms of only the low-energy 
observables, as we anticipated in eq. ([RJ). Such a ratio can be rewritten more explicitly: 



2 

I ~ / / / / -, \ ~ SIM ^f7.'W 

80, (37) 

l)ii,e \ m 2 / Sill Z(7 12 COS t7 23 



v m i)r/, _ /m 3 \ 2 sin20 23 



where for the estimate we have taken for 9yi and m 2 the values selected by the LMA solution. 
For the case of SMA solution that ratio increases to 10 3 — 10 4 . These results clearly hold 
in both scenarios (A) and (B). Therefore, this estimate can directly be translated into a 
prediction for the ratio of the decay rate of r — ► ^7 and \i — > as we shall show in the next 
Section. We also recall that since # 23 ~ 45°, the first generation is mixed with a state which 
is an equal (and indistinguishable) mixture of the flavour states v^v T . It is not surprising 
therefore that the ratio (nr- ) re / (nr~ )^ e be of order one: 

i-fp«tanfcs~l. (38) 

Notice also that the size of the lepton-flavour violating entries (m~)jj is about a factor 2 
larger in the scenario (B) due to the extra contribution driven by the Yukawa couplings Yz,s- 
This implies that the related decay-rates are further enhanced by a factor 4 in the scenario 
(B). Moreover, in the scenario (B), we have similar predictions for the sdown sector, namely: 

(m 2 - c ) 6s /m 3 \ 2 sin 2# 23 (m|) bd 
(mJ-J sd \m 2 J sin 20i2 cos 023 \ m %)sd 
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which show that the relative flavour-violation in the lepton and quark sector should be 
comparable in magnitude. 

Finally, we would like to comment also upon the case in which the 1-3 mixing is restored 
in the lepton mixing matrix. So for discussion we consider the present upper bound on 6 ri 
p8[ and take sin #13 = 0.1. This leads to an enhancement of the LFV entries (nx-)^ and 
(m|) re driven by the largest mass r/23, with respect to the #13 = 0, namely 



( m |; 


fie 




( m |; 


fie 


^13=0 ~ 


( m |; 


re 




( m |; 


re 


013=0 ~ 



\m 2 / sin 6^12 cos Byi 



1 - 



— —a a 10 - 40 

m 2 / Sill V\2 COS Vyi 



Therefore, since the entry (m|) TjU is not modified, the ratio of eq. ( |3?| ) becomes ~ 7, while 
that in eq. (|38|) remains the same. By considering the values for and #12 of the SMA 
solution, we find for sin #13 = 0.1 a much stronger relative enhancement, 

10 3 , 

-10 3 . (41) 



(m|; 


fie 


6»i3^0 


(m|; 


fie 


6»i3=0 ~ 


( m |; 


re 


013^0 


( m |; 


re 


6>i3=0 ~ 



Analogous results are obtained for the entries (m 2 r £ 



6 The £i — > ij^y decay rates 

Let us briefly recall here some points related to the computation of the ^ — > ^7 decay rate. 
The effective operator responsible for such a decay can be parameterised as 

£ eff = ig.ml&ll^i: + C^fan^F^ (42) 

where g e is the electromagnetic coupling and we use two-component spinor notation. This 
leads to the branching ratio 

BR{l L -> i n ) = (\&l\ 2 + |Cg'| a ) BR{t t -> l jVi v 3 ), (43) 

where in the specific cases the lepton-flavour conserving branching ratio are BR(fi — > 
ev^Ve) ~ 1, BR(t — > [iv T v^) ~ 17% and BRir — > eu T u e ) « 18%. For the numerical anal- 
ysis we have taken into account all contributions involving one-loop slepton-chargino and 
slepton-neutralino exchange, by using the complete formulas given for example in ref. |16| . 



Analogously to the case of the MSSM with singlets iV [IB], also in this scenario the main 
contributions come from tan /3-enhanced diagrams with chargino exchange. In the mass- 
insertion approximation, we recall that the parameter dependence of C % 1 (the dominant 
coefficient) is 

C l-T^— ^ tan/3, (44) 
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where m is an average soft mass. 

The main feature of our picture is the possibility to relate in a model-independent way 



the LFV of different sectors, as eqs. (|lq) and (|T7| ) demonstrate. Therefore, if we take the 



corresponding ratio of the BRs and take into account the estimate in fl3"T|), we find: 



BR(t — /i 7 ) / {ml) T A 2 BR(t -> ^i/ t i/ m ) 



10*. (45) 



For the case of SMA solution, this would become much larger, i.e. ~ 10 7 . Analogously we 
can predict 



BR(t -> erf) ( (m|) TC \ 2 5i?(r -> ez/ T z/ e ) 



-> e 7 ) \ (m| ) jUe / BR(ij, -> ez^z/ e ) 



10 _1 . (46) 



Let us now discuss the results obtained by a more detailed numerical study with exact 
solutions of the RGEs and exact diagonalization of matrices involved in the computation of 
the branching ratios. First in Fig. ||] we present the behaviour of the LFV parameters 5^ 

and Sfj defined as 

5 i = MM, 4 = MM, ( 47) 

L d c 

(m|,mj c are the average L and d c squared masses) as a function of the coupling A 2 in both 
scenarios (A) and (B), for two values of My, i.e. Mt = 10 11 GeV (upper left panel) and My = 
10 14 GeV (upper right panel) and for representative values of other parameters. For each 
scenario the upper and lower curve refers to and 5^ e (or S^ e , cfr. eq. fl38|)), respectively, 

and similarly for 5^ and 5f d , in the scenario (B). Due to the quadratic dependence in the 
RGEs on the LFV Yukawa matrices Yt,Ys,Yz, the S parameters scale as (A2)~ 2 . We 
notice that the size of LFV in the 2-3 and 1-2 sectors maintains a constant ratio, ~ 10 2 , 
independently of the scale (cfr. the estimate in eq. (|37D). In correspondence of each value 
of the scale Mt, there is a minimum value for A2, below which the RGE solutions blow up0? 
which is approximately A™ 111 ~ 3 x 10~ 4 (Mr/10 11 GeV). We can notice that in the scenario 
(B) Sjj gets larger than in case (A) by about a factor 2. Moreover, 5ff are about a factor 2 

smaller than 5^. This is due to the fact that in the evolution to low-energy the squark mass 
gets heavier than the slepton mass, because of the gluino- driving. This different increase 
is reflected mostly in the average scalar masses (which are determined at the SUSY scale) 
and to a less extent in the off-dagonal entries (m|)y, (m| c )jj (which instead are determined 
at the intermediate scale Mt)- The partial compensation of these effects explains why, for 
Mt = 10 11 GeV, for instance, the ratio 5 L /5 dc is smaller than 4, even though ~ 300 GeV 
and m d - c ~ 600 GeV. 

In the same figure we also show the behaviour of the branching ratios (lower panels) for 
tan/3 = 3, as a representative case. For each decay the lower and upper curve refers to 
the scenario (A) and (B), respectively. It is striking to notice that the constant-ratio rule 

12 One has also to check that A2 is not too large at Mt so that it remains perturbative up to Mq (cfr. 
the A 2 -RGE in eq. (Ell) in Appendix). 
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-i v f ? 0.2 0.1 0.6 0.8 1 1.2 1.1 



A 2 

Figure 2: In the upper panels we show the flavour-violation parameters S L and 5 dc as a function 
of the coupling constant A 2 (at M T ) for M T = 10 11 GeV (left) and M T = 10 14 GeV (right) in both 
scenarios (A) and (B). As for the neutrino parameters, we have taken #23 = 45°, 9u = 33°, #13 = 
and mi = 0, rri2 = 7 x eV, 7713 = 5 x 10~ 2 _eV. In correspondence of each scenario, the upper 
lines refer to 6^ and 5^ and the lower ones to 5^ e and 5g° d (the 5 dc are non-zero only in the case (B)). 

We recall that 5^ e = 5^ e as well as 8 d ° d = 8 d ° d . In the lower panels we show the resulting branching- 
ratios for the decays r — > ^7, fi — > e7 and r — > e7- The horizontal lines denote the experimental 
upper bounds of such BRs. For each BR the lower and upper curve is obtained in the scenario (A) 
and (B), respectively. In all panels, we have fixed = mo = 200 GeV at Mq] the corresponding 
average slepton mass at low energy is ~ 300 GeV (250 GeV) for My = 10 11 GeV (10 14 GeV). 
Finally, the parameters fixed at low energies are tan/? = 3, \i = 300 GeV and M2 = 180 GeV; the 
latter corresponds to M g = 490 GeV (290 GeV) at M G for M T = 10 11 GeV (10 14 GeV). We have 
also set Ai = A2 at Mt- 15 






< m d >/(1 TeV) 






tan/? = 3 




- "" M T =1 0" GcV ( 


vl 3 = 590 GcV) 



500 1000 1500 2000 2500 3000 3500 



m sol (GeV) 



Figure 3: In the upper panels we display the behaviour of the l{ — ► £,-7 BRs as a function of the 'left- 
handed' selectron mass in both scenarios (A) (left panel) and (B) (right panel) for My = 10 11 GeV 
and A2 = 1CP 3 (at Mp). For each BR the lower and upper line refers to the case with tan/3 = 3 
and 30, respectively. The horizontal dashed-line marks the BR(fj, — * cy) experimental bound. The 
gaugino mass and [i have been taken as in Fig. ^. In the lower panel we show the corresponding 
parameters 5fJ as emerge in the scenario (B), as a function of the selectron mass (solid and dashed 
lines). The dotted lines denotes the average mj c mass. The value of the gluino mass M3 as obtained 
at low energy is also indicated. 
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Figure 4: The same as in Fig. | for M T = 10 14 GeV and A 2 = 1. 
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displayed in eqs. (|45| - [Iq ) between the 5i? of the decays r — > ^7 or r — ► e7 with that of 
/i — > e7 is preserved in each case (A) and (B) as well as for any value of Mt, irrespectively of 
the value of A2. Moreover, in the scenario (B) with all the Yukawa couplings "Yt,s,z & t work 
these rates are larger since the LFV in the slepton masses is enhanced, as seen in the upper 
panels. The present bound on the BR(fi — * &y) (shown by dashed horizontal line) constrains 
A 2 to be larger than ~7x 1CT 4 and ~ 0.46 for M T = 10 11 GeV and 10 14 GeV, respectively 
in the case (A). In the picture (B) these bounds get a bit more stringent, A2 > 2 x 1CT 3 
and A2 > 0.65 for Mt = 10 11 GeV and 10 14 GeV, respectively. The T-induced seesaw strict 
predictions for the ratio of the BRs tell us that, in view of the future experimental sensitivity, 
both decays r — > ^7 and fi —>■ ey could be observed if the latter has at least a branching 
ratio of order 10~ 12 . For BR(fx — > ej) < 10~ 12 , we have that BR(t — > /ry) < 10~ 9 , below 
the expected sensitivity. On the other hand, if we depart from the limit #13 = and allow 
6>i3 as large as 0.1, we would get BR(r — > /ry) ~ 10 ■ BR(fi — > ej). As for BR{t — > erf) it is 
always predicted to be one order of magnitude smaller than BR(fi — > cy). 

In Fig. |3| we show BR(fi — > e-y), BR{r — > /ry) and BR{r — ► cy) as a function of the left- 
handed selectron mass at low-energy with Mt = 10 11 GeV and A2 = 10 -3 for both scenarios 
(A) (upper left panel) and (B) (upper right panel). For each BR, the upper and lower curves 
correspond to tan (3 = 30 and 3, respectively. In the scenario (A) for tan (3 = 3 there are no 
constraints on mg from the BR(fi — > erf) bound, while for tan /3 = 30 we have mg > 900 
GeV. More stringent lower bounds on m g can be deduced in the scenario (B). Therefore in 
the allowed range for mg, we have BR(r — > fij) < 2 x 10 -8 and BR(t — > ej) <2x 10 -12 . 
Again notice that the constant-ratio rule for the BRs is maintained also in this case where 
rrie is varied, confirming the fact that this rule depends only on the low-energy neutrino 
parameters and not on the details of the model, such_ as the soft-breaking parameters or 
Mt- In the lower panel, we show the behaviour of 5 dc and the average squark mass m 



(in- 



versus m g . The fact that S dc are about a factor 1.5 larger for tan/? = 3 as compared to the 
case with tan f3 = 30 is due to the combined increase for lower values of tan f3 of both Y u 
(see eq. |12])) and the top Yukawa coupling, Y t ~ mt/(v sin/3), which influences the running 



of Y v from low-energy up to Mt (see eq. (|5q ) in Appendix). We recall that the whole 
increase of Y u implies an increase of Yt as well of Y z and Y$ through the relation fl2~4|). 
Then the quantities Y^Y^ and YgY^ (scaling as 1/ sin 4 /5) trigger the flavour- violation in 
the matrix m d ~ c . The information on 5 dc and m d - c could be useful for the comparison with 
the present bounds on the flavour-violation parameters S dc extracted from the meson mixing 
measurements [ 29 1 . To this purpose we need to know also that the gluino mass at low-energy 
is M 3 ~ 600 GeV for M T = 10 11 GeV. For example, for tan (3 = 3 and m~ e > 300 GeV, we 
have rrijc ~ 600 GeV and 5f d { or 8^) ~ 4 x 10~ 4 which is below the limits from the K° — K° 
(or Bd — Bd) mixing parameter. 

In Fig. [| the same analysis has been performed for Mt = 10 14 GeV and we have chosen 
A2 = 1, in such a way that the ratio Mt/\i is the same as in the previous example. In this 
way, the size of the matrix Yy is the same at the scale Mt- All other parameters, such as M2 
at low energy, are the same as in Fig. |3|. Upon comparing with the previous case with lower 
Mt, we observe that for larger Mt the BRs are smaller by a factor 5 which is due to the 
smaller energy interval of the running, namely [log(Mo/10 n GeV)/ log(Mc/10 14 GeV)] 2 ~ 
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5. For the same reason, for Mt = 10 14 GeV the parameters 5 dc are smaller by a factor 2 or 
so. 



7 Conclusions 

The neutrino experimental observations pointing to sizeable lepton mixing have been en- 
couraging to further investigate the implications of lepton-flavour violation in extensions of 
the Standard Model. This work may be placed among these attempts. In particular, we 
have considered the SUSY seesaw mechanism obtained through the exchange of heavy SU(2) 
triplets. On comparing with the more popular seesaw scenario with the exchange of heavy 
singlets, our scenario is more predictive since the source of LFV at high energy, i.e. the 
Yukawa matrix Y^, can be directly connected to the low-energy observables, encoded in the 
coupling matrix Y u . The Yukawa Yt induces radiative corrections in the mass matrix of 
the sleptons L and as a result, even in the case of universal scalar masses at the GUT scale, 
LFV off-diagonal entries are generated. Therefore the flavour structure of the slepton mass 
matrix can be determined solely in terms of the low-energy neutrino parameters, i.e. the 
neutrino masses and mixing angles. The most remarkable feature of this scenario is that 
there is a rigid entanglement of the flavour-violation among different generations, as dis- 
played in eqs. (|3"T| - |3"8D , which does not depend on other details of the theoretical framework. 
This implies in particular that the ratio of the branching ratios of the decay r — ► /ry (or 
r — > ej) and /i — > cj can be predicted and turns out to be ~ 10 3 (or 10 _1 ). We have 



first derived these estimates by only taking into account the neutrino parameters (f4~5[-[46|) 
and have confirmed them by more detailed numerical computations, as shown in Figs. 2, 3 
and 4. Furthermore we have embedded this picture in a 'minimal' SU(5) scenario in which 
the triplet states T fill the 15-representation together with other coloured partners S, Z. In 
such a case flavour violation is also induced by radiative corrections on the mass matrix of 
the <i c -squarks. By imposing the GUT scale S , f/(5)-universality relation among the Yukawa 
couplings of those states to the matter multiplets, we find that the size of flavour violation in 
the lepton and quark sectors is comparable. This implies that, similarly to what happens in 
the lepton sector, the amount of flavour-violation between different quark sectors is strongly 
correlated. For example, one could predict the supersymmetric contribution to B s — B s 
mixing m terms of that to K° - K° mixing. It would be interesting to further explore this 
point and other implications of the SUSY T-induced seesaw. 
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Appendix 

In this appendix we first present the parametrization used for the chargino and neutralino 
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mass matrix since this determines the relative sign between the Yukawa and gauge terms 13 
in the renormalization group equations of the trilinear soft-breaking terms (see below (|55|)). 
Then we have determined the renormalization group equations in the MSSM with the 15 + 15 
representation of SU(5) at one-loop, which are therefore valid in the energy range between 
Mq and the triplet mass scale M T . 

The chargino mass matrix term is given by: 

-C ch = ( W+ \ T ( M 2 gv 2 W W~ \ +h.c, (48) 
\H+ J \ gvx fi ) \ H 2 J 

and that regarding the neutralino mass sector is: 



/ B \ 1 ( 
W° 
Hi 



Mi 




M 2 ^gv x 



7^ \ 



V2 



gvi 



o 



The renormalization group equations for the gaugino masses M a , (a 

,dM, 



( B \ 
W° 
HI 
V H« J 



+ h.c. (49) 



1, 2, 3) are: 



16tt 2 



dt 



Zg z a B a M a , 



(50) 



where the coefficients B a are given in eq. (|T9|) . The RGEs for the Yukawa couplings are: 



16tt z 



16tT 



16tT 



16tt z 



,dY T 
dt 

dYs 
dt 

dY z 
dt 

,dYe 

dt 
dt 



Y 7 



~gl - 7g 2 2 + YlY e + 6Y^Y r + Tr(Y^Y T ) + Y Z Y Z + |Ai| 



+(y^y: + y!y*)y t , 



Yc 



e e 

4 



g\ - \2g\ + 2Y*YJ + 8Y^Y 5 + Tr(Y+ Y s ) + 2Y Z Y T Z 



1 1 YtY e + 5Y^Y Z + Tr(Y^Y z ) + 3Y^Y T 



+2(Y d Y\ + Y Z Y\)Y S 
Yz 

L 1C 

J + 2Y S Y+ )Y Z , 

g\ - 3g 2 2 + 3(YtY e + Y^Y T + Y^Y Z + | A x | 2 ) + Tr(YtY e + 3Y+Y d ) 



1 2 _ 3 2 _ _16 ; 



+2(Y d Y d 

27 

15 
7 



Y d 



2 o 2 
r tv I Qvt 



16 



^ + 3(YtY d +|A 1 | 2 )+YtY M 



+Tr(YjY e + 3Y^Y d ) + 2(Y Z Y^ + 2Y 5 Y+ )Y„ 



13 We have indeed found some discrepancy in the literature. Our results, for example, are in agreement 
with the RGEs of |30) . In the latter work there is consistency between the sign of the gaugino mass terms 
and the RGEs of the trilinear terms provided a (missed) minus sign is accounted in front of the gaugino 
mass in the matrix of eq. (2.7) (or, equivalently, provided the i factor in the off-diagonal blocks is removed). 



20 



167T 

16tt 
16tt 



,dY u 

dt 
2 d\\ 



dt 

2 d\2 

~dT 



Ai 

A 2 



~9 2 i " H - fgl + 3YlY u + YjY, + 3|A 2 | 2 + 3Tr(Y^Y u ) 
7A 2 + Tr(Y T Y^ + 2Y e Yt + QY d Yl) - \gl - 7g. 
7\j + 6Tr(Y u Yt ) - ^gf - 7g 



(51) 



where Y^ and Y M are the Yukawa coupling matrix of the down and up quarks, respectively. 
Regarding the mass parameters of the superpotential, the RGEs are: 



2 dfi 

l07T — = a 

dt ^ 



3Tr(Yt Yu + YjY d ) + Tr(YtY e ) + S^) 2 + |A 2 | 2 ) - ^gf - 3g, 



16n 2 



,dM 7 



16tt 



16tt 2 



dt 
. dM s 

dt 
, dM z 

dt 



M T 
M s 
M 7j 



Tr(Y^Y T ) + | An_ | 2 + |A 2 | 2 - -g\ - 8g 2 2 

TrrvtV ^ 16 a 2 40 n 2 



1 



16 



Tr(Y^) - --g\ - 3g 2 2 - ^g\ 
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(52) 



The RGEs for the sfermion mass matrices are: 



16tt 2 



,dm 2 



dt 



m|(YlY e + 3Y^Y T + 3Y^Y Z ) + (Y|Y e + 3Y^Y T + 3Y^Y z )m| 
+2 (Y\m 2 H Y e + Ytm 2 - c Y e + 3Y^mf Y T + 3Y|m 2 - c Y z + 3Y^m|Y T 



+2>Y z m 2 z Y z ) + 2(A^A e + 3A^A T + 3A^A Z ) - -Mfgf - 6M 2 gl 

(2m| + Am 2 Hi ) Y e Yt + 4Y e m|Yt + 2Y e Yt m 2 - c + 4A e At - ^M^, 

2m 2 - c (Y,Yt + 2Y 5 Y^ + Y Z Y^ Z ) + 2(Y d Yt + 2Y 5 Y^ + Y z Y^)m 2 - c 
+4 (Y d m 2 Hl Yj + Y d m|Yj + 2Y 5 m 2 r rY 5 + Y z m|Y^ + 2Y 5 m 2 5 Y^ 
+Y z m 2 z Y z ) + 4(A d A* + 2A S A\ + A Z A Z ) - ^M 2 l9 j - ^M 2 gl 

(2m 2 - c + Am 2 H2 ) Y u Yt + 4Y u m|Yt + 2Y u Yt m 2 ~ c + AA U A[ 
-|M 2 , 2 - |m 3 2 , 3 2 , 

(m| + 2*4) Yt Y u + (m| + 2m 2 Hl ) Y^Y, + (yt Y u + YjY d ) m| 

+2Ytm 2 - c Y tt + 2Y\m%Y d + 2A\ A u + 2A d A d 

2 32 
--M 2 s 2 - 6M 2 g 2 2 - -M 2 g 2 3 . 
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,dm 2 



16tt 2 



dt 
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dm\ 
16tt 2 u 



16tt 2 



dt 



dm 2 



dt 
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The RGEs for the other soft-breaking masses are: 



_ „ 9 dm% 

16tt 2 — = 2 



m 2 ( | Ai| 2 + Tr(Y^Y T )) + 2Tr(Y^m|Y T ) + 2m 2 Hl \X 1 \ 
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+Tr(A^A T ) + |^i| 
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16tt 2 



, dm 

~dT 

,dm 2 s 

~dT 

16^ 
dt 

167r 2^l 

dt 



Mfg( - 16Mig z 2 , 
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2 (m||A 2 | 2 + 2m^|A 2 | 2 + |A 2 | 2 ) - ^-M\g\ - 16M 2 2 g 2 2 , 



2m|Tr(YT Y 5 ) + 4 [Tr(Y^ m%Y s ) + Tr(A f 5 A 5 ) 



§M 2 , 2 - ^M 3 2 , 3 2 , 



= 2m|Tr(Y^Y z ) + [Tr(Y^m 2 - c Y z ) + TY(Y* mf Y|) + Tr(A z A z ) 



16tt 2 



, <im 2 



2 32 

2 32 
--M^ 2 -6M 2 ^- y M^ 2 , 



16tt 2 



, dm 



Hi 



dt 



2m 



Tr(Y e Yt) + 3Tr(Y d Y*) + 6|Ai| 2 J + 2 ^(Yjm^Y,, + Y d m|Yj 



|2 2 



16tt 2 



Am 



h 2 



dt 



+Tr(Y e m|Yt + Yt m 2 - c Y e )] + 6^ 
+2 [3IV(A d Aj) + Tr(A e At) + 3|^!| 2 ] - ^Mfg 2 - QMjgl 

H 2 [3Tr(Y u Yt) + 6|A 2 | 2 ] + 6Tr(Yt m 2 - c Y u + Y tt m|Yt) 
+6|A 2 | 2 r4 + 6 [Tr(A u At) + |A 2 | 2 ] - ^M 2 g 2 - 6M 2 ^ 2 . 



(54) 



As for the soft-breaking trilinear coupling matrices we have: 
,dA T 



16vr 2 



dt 



= A 7 



YtY e + 9Y*Y T + 3Y Z Y Z + Tr(Y^Y T ) + |Ai| 2 - -g\ - 7g. 



+ (Y e T Y: + 9Y T Y^ + 3Y|Y* ) A T + 6A^Y|Y T 



+2Y, 



3Y^A Z + Tr(Y^A T ) + X\A 1 + \m i9 1 + 7M 2 g 2 2 

o 



167T 2 



dt 



= Ac 



12YJ Y 5 + 2Y*Yj + 2Y* Yf + ^(Yj Y s ) - - 12 a 



+2 (6Y 5 Y^ + Y.Y* + Y z Yi) A 5 + 4(A,Y^ + A Z Y^)Y 5 



+2Yc 



2Y*Aj + 2Y^Al + Tr(Y* A 5 ) + -M x g\ + 12M 2 ^ 2 



16tt 2 



,dA 2 
~dt 



A; 



Y|Y e + 3YyYr + 7Y^Y Z + Tr(Y^Y^) 



7 n n 16 2 

i - 3 #2 - y#s 



+2 (2Y 5 Y* + Y d Yj + 2Y Z Y^) A z + 6Y Z Y^A T 



16tt 



,dA e 



+2 
A P 



7 16 
AdY^ + 4A 5 Yj. + — M x g\ + 3M 2 g 2 2 + — M 3 ^ 3 2 



5YlY e + 3Y^Y Z + Tr(Y e Y| + 3Y d Y^) + 3|Ax| 2 + 3Y^Y T - \g\ - 3g. 



22 



,dA a 
~dt 



167T' 



16?r 



i dA u 
~dT 



,dA x 
~dT 



16tt 2 



,gL4 2 
~dT 



2YtA e + Tr(A e Yt + 3A d Y+) + 3A^ + 3Y^A r + 3Y^A Z 



+2Y e 

+ \m i9 \ + 3M 2 g 2 2 
5 



A n 



7 16 

5Y+Y, + YlY u + Tr(Y e Yt + 3Y d Y+) + 3|Ax| 2 - -g\ - 3g 2 2 - -g 2 3 



+2 (2Y S Y+ + Y z Yy A d + 4 (2A 5 Y^ + A z Yy Y d 
+2Y d [2Yj A d + Yt A u + Tr(A e Yl + 3A d Y+) + 3A^i 



7 16 

+— M^ x 2 + 3M 2 s 2 + yM 3 ^ 2 



A. 



13 



16 



5Yt Y u + Y+Y d + 3Tr(Y„Yt) + 3|A 2 | 2 - —g{ - 3g 2 2 - -g\ 



15' 



+2Y U 



13 



2Yt A u + Y\ A d + 3Tr(A„Yt ) + 3A*A 2 + —M x g\ + 3M 2 ^ 2 



16 



15 



M 3 g 2 3 



A, 

+2A :l 
A 2 

+2A 2 



2Tr(Y e tY e + 3Y+Y,) + Tr(Y^Y T ) + 21|A!| 2 - -g\ - 7g 2 2 



Tr(Y^Ar) + 2Tr(YtA e + 3Y+A,) + \M x g\ + 7M 2 g 2 2 

5 



6Tr(YtY u ) + 21|A 2 | 2 -^ 2 -7^ ; 

5 

6Tr(Yt A u ) + \m i9 2 x + 7M 2 <? 2 
5 



(55) 



Clearly, all the RGEs shown above are valid in both scenarios (A) and (B). In the former 
scenario, the condition Y$ = Y z = at Mq ensures that Yukawa couplings Y s, Y z and the 
parameters As, Az are not radiatively induced, therefore they can be simply switched off in 
the r.h.s of any RGEs. Beneath the mass scale Mt, we recover the RGEs of the MSSM by 
switching off Yt,s,z, Ai, 2 , m TT , m 2 s m 2 z ^ and Ax t s,z, For completeness, we report also 
the RGE of the d = 5 neutrino-operator Yukawa matrix Y u valid below Mt in the MSSM 



167T 



,dY u 
dt 



6 



g\ - Qg 2 2 + 6Tr(Yt Y M ) + Y„Y^Y e + (Y^Y e ) T Y y . 



(56) 
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